Accepted by The Astrophysical Journal 

Preprint typeset using I^TeX style emulateapj v. 11/10/09 



EVOLUTION OF BLISTER-TYPE HII REGIONS IN A MAGNETIZED MEDIUM 

Leo Gendelev^ and Mark R. Krumholz ^ 

Department of Astronomy, University of California, Santa Cruz, CA 95064 
Accepted by The Astrophysical Journal 

ABSTRACT 

We use the three-dimensional Athena ionizing radiation-magnetohydrodynamics (IRMHD) code to 
simulate blister- type HII regions driven by stars on the edge of magnetized gas clouds. We compare 
these to simulations of spherical HII regions where the star is embedded deep within a cloud, and to 
non-magnetized simulations of both types, in order to compare their ability to drive turbulence and 
influence star formation. We find that magnetized blister HII regions can be very efficient at injecting 
energy into clouds. This is partly a magnetic effect: the magnetic energy added to a cloud by an HII 
region is comparable to or larger than the kinetic energy, and magnetic fields can also help collimate 
the ejected gas, increasing its energy yield. As a result of these effects, a blister HII region expanding 
into a cloud with a magnetic field perpendicular to its edge injects twice as much energy by 5 Myr as 
a non-magnetized blister HII region driven by a star of the same luminosity. Blister HII regions are 
also more efficient at injecting kinetic energy than spherical HII regions, due to the recoil provided by 
escaping gas, but not by as much as predicted by some analytic approximations. 
Subject headings: HII regions; ISM: clouds; ISM: kinematics and dynamics; ISM: magnetic fields; 
MHD; radiative transfer; stars: formation 



1. INTRODUCTION 

In recent years it has become evident that HII re- 
gions - formed when massive stars ionize their surround- 
ings - are of crucial importance to the evolution of gi- 
ant molecular clouds (GMCs) since they photoevaporate 
their nascen t clouds, tri^^er star formation, and drive 
turbulence. |McKee fc Williams ( 1997| ) find that mas- 
sive stars inject much more energy into the interstel- 
lar medium (ISM) by emitting ionizing radiation than 
through supernovae. Furthermore, semi- analytic mod- 
els show that HII regions can disrupt GMCs entirely, 
through a combination of photoevaporation and mechan- 
ical disruption, over a period of a few cloud crossing times 
- roughly 20 - 40 Myr ( Williams fc McKeeJ997; Matzner, 
2002[ iKrumholz et al.||2UU6[ iGoldbaum et al.||:^(jll| . As' 
the ionization front moves tnrough the cloud, it sweeps 
up neutral gas, potentially increasing the star formation 



rate in the dense shell (Elmegreen & Lada" fT977liJWhit- 
[worth et al. ( 1994 ) ; however see [ Dale et al. ( 20(J7a ), who 
argues that this effect is generally small.) HII regions 
may also be important for driving turbulence in GMCs. 
Observed CMC lifetimes are significantly longer than 
the free-fall t ime scale ass ociated w ith the d issipation 
of turbulence (IK rumholz fc Tan 2007l|Fukui et al. 2009). 
|Matzner| ( 2002 ) finds that HII regions provide more en- 
ergy for turbulence than th e combined effects of st ellar 
winds and supernovae, and Krumholz et al. (2006) and 



Goldbaum et al. (2011) argue that this energy is suffi- 
cient to drive turbulence over observed CMC lifetimes. 
Hence it is important to study HII regions and their en- 
ergy injection mechanism in greater detail. 

Recent numerical studies have examined HII regions in 
turbulent media as a possible mechanism to explain the 
qualitative f eatures and star formation rates observe d in 
GMCs (e.g . Mellema et al.| l2006| [Dale et al.|[2Q07b| and 
Gritschneder et al. 2009). I'hese papers show that the 



bulent gas has important effects. In Gritschneder et al. 
( 2009 ) the turbulent energy is significantly increased (up 
to a factor of 4 in the cold gas) as the ionization heats 
the gas along channels of low density, compressing higher 
density gas into filaments with gravitational collapse oc- 
c uring in the tips of the pillar like structures. 



interaction of~the HII region with the pre-existing tur- 



Krumholz et al. (2007) perform the first ever numerical 
study of the expansion of an HII region into a magnetized 
gas. It is important to study the effects of magnetic fields 
since the magnetic energy in GMCs is comparable to the 
kinetic and gravitational energies (Crutcher|1999l). Over 
time the ionization front slows down enough due to the 
resistance by the magnetic field lines in the perpendicu- 
lar direction that the fast magnetosonic wave outruns it, 
disturbing neutral gas ahead of the ionization front, so 
that there is swept up material in between the fast mag- 
netosonic wave and the ionization front. In this region 
the stretched magnetic field acts as an energy reservoir. 

In this paper we expand on this work by examining 
magnetized blister- type HII regions, also known as cham- 
pagne fiows, which form when an ionizing star is situated 
towards the edge of th e CM C. A sketch adopted from 
Krumholz & Matzner] ( |2009) is shown in Fig. [T| com- 
paring blister and synimetnc HII regions. In the olister 
case the star is situated next to the edge of the CMC, so 
the ionization front will eventually reach the edge of the 
cloud and burst a hole through which hot ionized gas 
will be able to stream out at supersonic velocities into 
the low-density ISM. The b lister scenario was first envi- 
sioned by [Whitworth (1978) as a mechani sm to disperse 
molecular gas clouds and was studied by | Tenor io-Tagle 
( |1978l ) to explain some features of observed nebulae. In 
contrast with the symmetric case, the ionized gas within 
the HII region is not confined to the HII region, leading 
to an increased expansion rate and kinetic energy of th e 
ionization front due to the "rocket effect" ( Kahn p"954 
In the symmetric case the gas is confined within the H. 
region and cannot rocket away, so this effect does not ap 
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Figure 1. A sketch comparing th e bhster and symmetric HII 
regions ( [Krumholz &: Matzner|2009| >. 

ply. Hence we would expect more energy to be injected 
into the GMC by the blister HII region. It might seem 
coincidential that a star would be positioned close to the 
edge of a molecular gas cloud, but observational evidence 
points to th is being the rule rather than the exception 
(Israel 1978). This makes sense given that GMCs are 
turbulent - the turbulence creates a filamentary struc- 
ture in the GMC so that any star has a high probability 
of being born near an "edge" . There are only a few ex- 
amples of blister HII regions for which the mag netic field 
can be measured. One is Ml 7, in which Pelleg rini et al. 
( 2007 ) find that the magnetic field is strong enough that 
it has es sentially h alted the expansion of the HII region. 

Since Krumholz et al.| (|2007j) there have been sev- 
eral n umerical studies to take MHD effects into ac- 
count. [Henney et aD ( |2009D and iMackey fc Lim| ( [2010 ) 
perfornied IKMHD (ionizmg radiation magnetohydro- 
dynamical) simulations of magnetized globules. Given 
strong initial magnetic fields, |Henney et al.| (j2009) ob- 
served substantial deviations from symmetry, m partic- 
ular when the initial magnetc field was oriented perpen- 
dicular to the direction of ionizing radiation. The ionized 
sras as a whole was found to be dominated by magnetic 
pressure, and evidence was found that magnetic effects 
might be important in the formation of bright, bar like 
emission features in H II regions, a result confirmed by 
Mackey fc Lim (2010). Even more recently Arthur et al. 
( |201l| ) perfornied siniulations of HII regions expanding 
into a magnetized gas with turbulence, resulting in a 
morphology of striking similarity to observed HII regions 
with clearly destinguishable pillar like structures. How- 
ever, none of these studies have dealt with blister HII 
regions, and it would be difficult for those of them that 
focus on turbulent regions to do so, because these sim- 
ulations use periodic boundary conditions which cannot 
account for gas being "blown off" the computational grid 
- and in a turbulent medium there is always the oppor- 
tunity for material to be blown off. These studies are 
useful in that they provide evidence for the importance 
of magnetic fields and the ability of turbulence to re- 
produce observed features of HII regions, but they need 
to be followed up by the examination of the energetics 
rather than the qualitative consequences of MHD turbu- 
lence and a comparison of the energy injection efficiency 
of the symmetric and blister-type HII regions. 

Our work is the first numerical study of blister type 
HII regions evolving in the presence of uniform magnetic 
fields but without turbulence and is therefore comple- 
mentary to previous work on symmetric HII regions in 



a turbulent medium with and without MHD. The struc- 
ture for the rest of the paper is as follows: in ^ we 
present our computational approach and parameters, in 
^we present our results from the various types of runs, 
looking at symmetric and blister-type HII regions both 
with magnetic fields of varying orientations and without 
magnetic fields, in Q we discuss the results and com- 
pare to analytic approximations, and in ^we draw our 
conclusions. 

2. COMPUTATIONAL METHODOLOGY 

We perf orm the simula tions with the Athena 3D MHD 
code (Stone et al."2008'). Athena is a grid-based, static 
mesh refinement (SMR) code designed for simulating var- 
ious astrophysical MHD processes. It uses higher order 
Godunov methods, which are particularly efficient when 
used with static or adaptive mesh refinement, in combi- 
nation with the constrained transport technique, which is 
used to ensure that the magnetic divergence is preserved 
to machine precision. In additi on, we employ the radia - 
tion scheme first introduced in 



We run all simulations on the Pleiades cluster ai 
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(2007). 



at 256^ resolution using 64 processors for an average wall- 
time of 3 days per simulation. 

2.1. The Ideal MHD Equations 

Athena solves the equations of ideal MHD. In conser- 
vative form, they are: 
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where p is the density and v is the velocity of the gas, 
B is the magnetic field, P* = P + • B/2) is the total 
pressure, P is the gas thermal pressure, E is the total 
energy density, and pn is the density of the neutral gas. 
Equation ([T]) is the continuity equation (conservation of 
mass), ([2| is the conservation of momentum, where we 
have usea the approximation that our fluid has no viscos- 
ity (i^ = 0), (|3| is Faraday's law, or the induction equa- 
tion, where we have set the magnetic diffusivity to zero 
(t^ = 0), (|4]) is the conservation of energy equation, where 
Q and C are the radiative heating and cooling terms, re- 
spectively, ([5| is the continuity equation for neutral gas, 
and says that a change in the mass of neutral gas can 
only come through advection into other cells or through 
recombinations (7^) and ionizations (X), and (|6| comes 
from the non-existence of magnetic monopoles. 

2.2. Heating, Cooling, and Ionization Chemistry 

All heating, cooling, and ionization terms are as in 
Krumholz et al. (2007), where the heating a nd coo ling 



rates are adopted from. Koyama fc Inutsuka] ( 2002 ) for 
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neutral gas, and from Osterbrock ( 1999 ) for partially ion- 
ized gas. The heating and cooling equations are provided 
below: 



g = e^auH ^ 



47r|x - 



-(^'^-)+n^rK/, (7) 



C = ArecnenH+ + Aff{T)nenH+ + AKi{T)njj 

^Aiine{T)nenH+, (8) 

where n//, n//+, and rig are the number densities of neu- 
tral hydrogen, ionized hydrogen, and electrons, respec- 
tively, is the thermal energy added per ionization, a 
is the cross section for absorption of a photon at the ion- 
ization threshold for neutral hydrogen, Sn and are the 
ionizing photon luminosity and position of the nth star, 
X is the position where the ionization is taking place, r 
is the optical depth to ionizing photons between x, and 
Xn, given by 



r(x,Xn) = {(TnH)dl 



X 



(9) 



A^ec, A//, and Aune are the cooling rates due to recombi- 
nation radiation, free-free emission, and metal line emis- 
sion in the ionized gas, and Tki and A^i are the neutral 
ISM heating a nd cooling rates, computed using the ap- 
proximation of Koyama & Inutsuka (2002). 

We should point out a particular feature of heating and 
cooling rates that apply in the neutral medium (where 
= 0)- these particular heating and cooling curves 
have the property that there is a two phase equilibrium. 
At a given pressure there are two solutions: a high tem- 
perature - low density solution, and a low temperature - 
high density solution. We use this property of the neu- 
tral gas to initially set up our blister type problems, by 
placing the high density and low density halves of the 
computational domain in pressure equilbrium. 

2.3. Problem Setup 

We set up a rectangular grid which runs from -25.0 to 
25.0 pc in all 3 directions. The computational domain 
has outflow boundary conditions. In all runs we place a 
star at the center of the grid, with (x, z) = (0, 0, 0). We 
then use two possible initial density distributions. In the 
spherical case, we initialize the grid to a uniform initial 
density n = 63.0 H atoms cm~^ across the whole domain. 
For blister-type runs, we split the grid into two halves: all 
cells with X < are of higher density, nieft = 63.0 cm~^, 
while all cells with x > are of a low density, nright = 
0.055 cm~^. Given our cooling curves, the equilibrium 
temperature at a density of n = 63.0 cm~^ is T = 55.0 
K (sound speed cq = 5.74 x 10^ cm s~^), while at a 
density of n = 0.055 cm"^ it is T = 6.3 x 10^ K (sound 
speed Co = 6.14 x 10^ cm s~^). Thus the two sides are in 
pressure balance. The sound speed inside the HII region 
is Q = 8 X 10^. The mean particle mass for the neutral 
gas is 2.3 X 10~^^ grams. For the spherical MHD and the 
blister MHD runs we thread the initial magnetic field 
through the domain in the x direction. In addition, we 
perform one blister-type run with the initial magnetic 
field at 45° to the x axis and finally one run with the 



magnetic field in the y direction. In all MHD runs, the 
initial magnetic field is = 3.0 x 10~^ in the code units, 
which differ from cgs units by a factor of a/Stt, so that in 
cgs it is 10. 6/iG. A summary of some of these parameters 
is provided in Tab. [T| 

We parameterize the luminosi ty of our centra l source 
in terms of its Stromgren radius ( Stromgren |1939i , which 
is defined as: 
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where s is the ionizing luminosity of the star, jj^h and p 
are the mean mass per hydrogen atom and the density 
of the gas, respectively, and a^^^ is the recombination 
coefficient. We set = 1.5 pc for al the simulations, 
computed using the value of p for the dense half of the 
grid in the blister cases. The corresponding ionizing lu- 
minosity is 5 = 5.3 X 10^^ s~^, appropriate for a star of 
spectral type BO. 5. 

We configure Athena with the Roe solver based on the 
Godunov scheme in conjunction with the ctu integra- 
tor to produce the most accurate results, and we enable 
h-correction in order to eliminate carbuncle problems. 
We resort to using first order fiuxes because higher or- 
der fiuxes proved unstable at 256^ resolution in at least 
some of our runs. All other parameters relevent to our 
problem required to configur e Athena so as to reproduce 
our results are in Kru mholz et al.j ( j2007t ). 

A note on the blister type HII region setup: in reality it 
would clearly be a coincidence if the star was positioned 
right on the edge of the cloud. Blister-type HII regions 
are likely to form when a star is close to, but not directly 
on, the edge of the cloud. However ours is an instructive 
limiting case, since it is much harder to interpret and 
make sense of computational data that would result from 
a more realistic setup. It would be useful to extract as 
much information as possible from our idealized setup, 
and in a later paper expand our investigation to compare 
to the mor e realistic scenario such as Ml 7 (Pellegrini| 
irir]|2QQ7| )). 



3. RESULTS 

3.1. Symmetric Simulations 

We begin our analysis by revisiting the extensively 
studied classic - the symmetric HII region - where the 
star is situated deep within the cloud. 

3.1.1. Hydro Run 

As expected, in the absence of a magnetic field, the 
ionization front expands in a spherically symmetric shell 
as shown in Fig. [2] We can see that the density of the 
shell increases over time and that virtually all the kinetic 
energy is contained within the thin shell that bounds the 
ionization front. 

3.1.2. MHD Run 

In the presence of the magnetic field, the expansion of 
the ionization front is strongly suppressed perpendicular 
to the magnetic field lines, so that over time the symmet- 
ric HII region assumes the shape of a football, as seen in 
Fig. |3] From row 3, it is apparent that at 0.5 Myr - 
corresponding to the first column of the figure - the fast 
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Table 1 

Problem Setup 



Simulation Name 


Type 


B-Field Setup 


B-Field 


nieft 






Tleft 


^right 




Hydro 


Symmetric 


None 


None 


63.0 cm" 


-3 


63.0 cm-3 


55.0 K 


55.0 K 




MHD 


Symmetric 


X direction 


10.6 iiG 


63.0 cm" 


-3 


63.0 cm-3 


55.0 K 


55.0 K 




Blister-hydro 


Blister 


None 


None 


63.0 cm" 


-3 


0.055 cm-3 


55.0 K 


6.3 xlO^ 


K 


Blister-mhd 


Blister 


X direction 


10.6 iiG 


63.0 cm" 


-3 


0.055 cm-3 


55.0 K 


6.3 xl03 K 


Blister-mhd-vert 


Blister 


y direction 


10.6 ^iG 


63.0 cm" 


-3 


0.055 cm-3 


55.0 K 


6.3 xlO^ 


K 


Blister-mhd-45 


Blister 


45° 


10.6 ^jlG 


63.0 cm" 


-3 


0.055 cm-3 


55.0 K 


6.3 xlO^ 


K 
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Figure 2. Slices in the z = plane taken from the hydro run. The 1**, 2^^, and 3^^ columns correspond to 0.5, 1.5, and 5.0 Myr into 
the simulation, respectively. The 1"^* and 2^^ rows display the density and kinetic energy density of the hydro run, respectively. 



magnetosonic wave is just beginning to take the lead in 
front of the sheh. By 1.5 Myr this is one of the most 
prominent features of the plot. Most of the added mag- 
netic energy is contained in this region, implying that 
the magnetic energy injected into the cloud by the HII 
region increases at late times. 

3.2. Blister-Type Simulations 
3.2.1. Blister-hydro Run 

First we look at the blister-type scenario where there 
is no magnetic field present, in order to be able to better 
understand what effects the addition of a magnetic field 
has on the HII region. 

Density and kinetic energy slices of the computational 
grid are presented for 0.5, 1.5, and 5 Myr of this run in 
Fig. [4] Initially the expansion into the dense half resem- 
bles the expansion in the symmetric case - the ioniza- 



tion front shell is almost identical to the left hemisphere 
of the spherical shell in the symmetric non-mhd simula- 
tion (Fig. [2]). Over time, however, the deviation from 
symmetry becomes increasingly apparent. By 1.5 Myr 
there are slivers of the dense shell that extend further in 
the y direction. This effect is easiest to see from the ki- 
netic energy plot in row 2 of the same figure. Unlike the 
embedded case, there is a jet of gas blowing out of the 
cloud. Although this low-density material covers a wider 
area than the dense shell, its average kinetic energy den- 
sity is orders of magnitude lower than the kinetic energy 
density within the shell, so its kinetic energy is virtually 
negligible. 

3.2.2. Blister-mhd Run 

We now present the blister-mhd results that are the fo- 
cus of this paper. The magnetic field orientation we look 



Evolution of Blister- Type HII Regions in a Magnetized Medium 



5 



20 1 



-21 



10 i 



u 

Q_ 



! 



-10 ! 

-20 I 
20 



-22 
-23 

-24 

-25 
-10 



CP 

o 



10 



a 

CL 



O 



( ) 



I 

u 

-11 ^ 



-10 1 



-20 



20 I 



10 i 



u 

Q_ 



-10 1 



-20 i 



J 
1 



-20 -10 10 20 
r (pc) 



-20 -10 
r (pc) 



20 -20 - 



10 
r (pc) 



10 



20 



O 



-12 
-1 1 



o 



-12 



1-13 



Figure 3. The mhd run - same as Fig. p] except for the addition of the third row, which shows the change of the magnetic energy over 
time, AEb, where AEb = ((B^)^ + (By)^+ (Bz)'^)/2 - ((B^q)^ + (^yo)^ + (^^p)^)/2, where B^^,By^,B^^ are the initial magnetic field 
strengths in the x, y, and z directions, respectively. Note that within the HII region interior AEb < 0? which is why we plot the log of its 
absolute value. Everywhere else in the computational domain AEb > 0. The white lines in the top row represent the magnetic field lines, 
drawn from footpoints equally spaced along the left side of the region shown. 

in the blister- hydro case (row 2 of Fig[5|. However, the 
shell is more oblate, and in contrast to the blister- hydro 
shell, the kinetic energy decreases much more slowly with 
time. Thus one of the most important MHD effects in 
the blister-mhd case is that the magnetic field changes 
the nature of the expansion over time (as we will see in 
^3.3), by collimating the jet of gas streaming out of the 
cloud. 



at first is particularly useful since it is easy to compare 
and contrast with the blister-hydro run. The expansion 
of gas - both in the dense and low density portions of the 
computational domain - is suppressed in the directions 
perpendicular to the magnetic field (Fig. |5|, and the 
magnetic field streamlines the ionized gas blowing out of 
the cloud. 

A chunk of hot gas resembling a bullet bursts out of 
the cloud and by 1.5 Myr has already reached the edge of 
the grid, implying a speed of 15-20 km/s, so a significant 
amount of kinetic energy is both gained and lost over the 
course of the simulation. 

The kinetic energy in the dense half is all concentrated 
in the shell, which is very similar in structure to the shell 



3.2.3. Blister-mhd-vert Run 

Initially, similar to the blister-mhd case, a jet of gas 
bursts out of the cloud as seen in Fig. [6] However the jet 
is not streamlined and is not propelled at high velocities 
through the low-density medium. Its motion is highly 
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Figure 4. The blister-hydro run: the 1** row are the density shces while the 2^^ row contains the kinetic energy density. The density is 
plotted on a log scale so that effects in the low-density part of the computational domain can be distinguished. 



suppressed by the magnetic field perpendicular to it. The 
magnetic field gets stretched the most out of any of the 
runs. 

The kinetic energy of the cloud is concentrated in the 
shell just as for the other runs, but the magnetic field's 
suppression of the shell's motion in all but the vertical 
direction transfers most of the kinetic energy from the 
spherical part of the shell to the slivers of dense gas pro- 
truding into the low density medium by the end of the 
run {2^^ row). The magnetic energy density is plotted in 
row 3 of Fig. [6j As expected, since the ionization front 
is slow to expand perpendicular to the magnetic field 
lines, the fast magnetosonic wave disconnects from the 
front earlier than in any other simulation, resulting in a 
very prominent density fiuctuation leading the ionization 
front. 

3.2.4. Blister-mhd-45 Run 

This simulation is clearly a blend of the properties of 
the blister-mhd and blister-vert runs. As seen in Fig. [7| 
the magnetic field lines limit the expansion of the front 
perpendicular to them, but do not prevent the front from 
travelling at a sizeable velocity parallel to them. The 
result is an HII region that has rectangular structure. 

3.3. Comparison of the Symmetric and Blister 
Simulations 

In this section we compare properties such as the ra- 
dius, kinetic, magnetic, and total energies of the various 
simulations. 



3.3.1. Shell Radius and Mass 

The radius of the shell as a function of time can be 
estimated analytically using conservation of momentum 
and some simplifying assumptions. If we assume that the 
pressure inside the HII region is vastly dominant over the 
ambient pressure in the neutral gas into which it is ex- 
panding at all times, and that the density within the HII 
region is approximately uniform, then we can obtain an 
equation of motion for the shell from momentum con- 
servation. This equation of motion can then be solved 
by using a power-law similarity solution. We provide a 
derivation for both the s ymme tric and blister cases in 
Appendix |A| In S pitzer] ( 1978 ) there is a derivation of 
the symmetric case performed m a slightly different way, 
yielding a very similar solution. The solutions we find 
for the two cases differ only by a factor of 2^/'^, and are 
provided below: 



It 



I2t. 



4/7 



(spherical). 



and 



(- 



It 



^sh = Ts I ) (blister), 



4/7 



(11) 



(12) 



where ts = rs/ca. 

How valid are the approximations we used to derive 
the shell radius (Eqn.'s ( pT| and (12))? In reality, even 
though the assumption that the arnment pressure is neg- 
ligible in comparison with the pressure inside the HII 
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Figure 5. The blister-mhd run. All panels are the same as in Fig. [3] 

region is quite accurate at early times, the HII region in- 
ternal pressure decreases as the HII region expands, and 
so it becomes less accurate. Furthermore, even though 
the assumption that the density is uniform within the 
HII region for the symmetric case is valid, it is not nec- 
essarily quite as true for the blister case where gas is 
free to stream out of the HII region into the interstel- 
lar medium. Hence the radius as a function of time will 
likely not match the analytic solution perfectly for all 
times, especially for the blister case. This is discussed in 
more detail in Q 

We define the shell radius as the average radius of all 
cells whose density is greater than or equal to 1.0 Ipo 
(where po is the initial density of the neutral gas) in or- 
der to avoid taking into account the neutral, undisturbed 
gas and the low density medium, as well as to take into 



account the contribution to the radius from the gas in 
the fast magnetosonic wave which can have a density of 
Po < p < l.lpo- If we were to use a cut-off of say, l.lpo, it 
would make little difference for the hydro runs, since they 
are characterized by very thin, high-density shells. For 
the MHD, runs, though, using a cut-off of l.lpo would 
not only produce a smaller radius, it would lead us to 
severely underestimate the added magnetic energy, since 
most of this energy is contained in the mildly overdense 
region between the ionization front and the fast magne- 
tosonic wave. 

The shell expands slower for the hydro case than the 
analytic approximation (Eqn. (pT|). The average ra- 
dius in the mhd run is greater tHan in the hydro run 
due to the fast magnetosonic wave leading the ioniza- 
tion front (Fig. Isl). The radius of the shell is defined 
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Figure 6. The blister- vert run. All panels are the same as in Fig. [s] 

using the same criteria as in the hydro run, but this 
time the definition encompasses the material contained in 
between the ionization front and the fast magnetosonic 
wave since this material has a higher density than the 
neutral background. As expected, all the blister-type 
runs expand faster than their symmetric counterparts. 
Although there is some variation between the magnetized 
blister runs in the early-mid stages of the simulations, at 
late times they all converge with one another. This is 
due to the fact that, at late times, the radius is effec- 
tively set by the fast magnetosonic fiuctuation leading 
the ionization front. Thus in all the MHD cases, the ra- 
dius is determined not only by the speed of the actual 
ionization front, which varies between the cases, but also 
by the fast magnetosonic speed, which does not. 
We also calculate how much mass is swept up in the 



dense shell (as well as in the fast magnetosonic wave) and 
present these results in Fig. [9) The dashed lines repre- 
sent the hydro and mhd runs divided in half. Although 
the symmetric runs sweep up more mass than their blis- 
ter counterparts, the difference is less than a factor of 
two, so the fact that the blister runs send hemispheri- 
cal rather than a spherical shell into the cloud is partly 
compensated for by the fact that the shell expands more 
rapidly. The blister-mhd run sweeps up nearly twice as 
much mass by 5 Myr as the blister-hydro run, so the 
magnetic fields make a large difference with respect to 
sweeping up of mass. The orientation of the initial mag- 
netic field makes little difference in terms of sweeping up 
of mass for the blister-type runs. 

3.3.2. Kinetic Energy 
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Figure 7. The blister-mhd-45 run. All the panels are the same as in Fig. [s] 



We plot the kinetic energies of the various runs in Fig. 
TQ| and the specific kinetic energy, defined as the kinetic 
energy divided by the swept up mass, in Fig. pTl Note 



that, since we only include cells with p > l.Olpo in this 
computation, we are only evaluating the kinetic energy 
imparted to the dense cloud, and not the kinetic energy 
carried by the outflowing gas or deposited in the low 
density medium. 

The result that the blister-hydro kinetic energy is lower 
than the hydro kinetic energy is contrary to what we 
would expect to see based on the analytic solution. From 
Eqn's. and (12), we see that the predicted blister 



radius is larger by a factor of 2^/^ than the symmetric 
radius for any fixed time. The kinetic energy of the shell 
is given by Eke = l/2Mr^^, where M is the mass of 
the shell. The mass of the shell increases as r^^, so Eke 



should be larger by a factor of (1/2) (2^/^) (22/^)2 = 2^/^ 
for the blister- hydro run, where the 1/2 term is included 
to account for the fact that the blister-hydro shell is a 
hemisphere rather than a sphere. However, the spheri- 
cal part of the blister-hydro shell actually expands sig- 
nificantly slower than the analytic solution predicts (see 
§4|, so if most of the kinetic energy is concentrated in 
tEe spherical part of the shell, the total kinetic energy 
for the blister-hydro run could be lower than for the hy- 
dro run. This is indeed the case. As seen in Fig. |4]row 2, 
the slivers of the shell that expand in the y direction have 
very little kinetic energy density. Thus, even though they 
contribute significantly to the radius of the blister-hydro 
shell, the total kinetic energy is lower in the blister-hydro 
case than in the hydro case. 
In general, the non-mhd runs have both more total and 
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Figure 8. The mass averaged radii of all the runs. The blue- 
dashed curves represent the analytic solutions for the blister (upper 
curve) and spherical (lower curve) cases respectively. 
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Figure 9. The mass of the various runs summed over cells whose 
density is greater than 1.01 times the initial density. The dashed 
curves represent the hydro and mhd curves divided in half for com- 
parison. 

specific l^inetic energy than their magnetic counterparts. 
The total and specific kinetic energy is very sensitive to 
the orientation of the initial magnetic field. The blister- 
mhd run has more total and specific kinetic energy than 
the mhd run, but the mhd run has more total kinetic en- 
ergy and a comparable amount of specific kinetic energy 
to the blister- mhd- vert and blister- mhd-45 runs. It is in- 
teresting to compare the effect of going from the mhd to 
the blister-mhd run to that of going from the hydro to 
the blister-hydro run. Going from symmetric to blister 
has the effect of increasing both the total and specific 
kinetic energies in the presence of a magnetic field, but 
has the opposite effect in their absence. It is also inter- 
esting to note that the blister-hydro run has less specific 
kinetic energy than the hydro run. This is not surpris- 
ing, since the mass swept up by the blister-hydro run is 
significantly larger than half of the mass swept up in the 
hydro case (Fig. [9|. 

3.3.3. Magnetic Energy 



Figure 10. The total kinetic energy of all the runs. The kinetic 
energy is calculated only considering cells whose density is greater 
than 1.01 times the initial density. The dashed curves represent 
the hydro and mhd curves divided in half. 
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Figure 11. The specific kinetic energy of all the runs calculated 
only considering cells whose density is greater than 1.01 times the 
initial density. 

For the magnetic energy, we expect the results to be 
opposite of those for the kinetic energy. The runs with 
the least kinetic energy should actually have the most 
magnetic energy since energy that does not go into mo- 
tion is instead stored as distortions of the magnetic field. 
We present these results in Fig. [12] The dashed line is 
the mhd run divided in half for comparison. 

We calculate the change in magnetic energy using the 
same criterion as for the kinetic energy (p > l.Olpo) 
not only to avoid taking into account fiuctuations that 
occur outside of the cloud, but also to exclude cells on 
the edge of the grid. In our simulations the fast mag- 
netosonic wave eventually reaches the edge of the grid, 
stretching the magnetic field lines outside of the compu- 
tational domain. If we were to take into account edge 
cells, we could see a net loss of magnetic energy as the 



1x10"*^ 



8x10^ 



6x10^ 



4x10^ 



2x10^ 



I I ' ' ' 

• •MHD 

• •MHD- HALF 

• • BLISTER -MHO 

• • BLISTER-MHO-VERT 

BLlSTER-MHO-45 



Evolution of Blister-Type HII Regions in a Magnetized Medium 

10 



11 




2 3 
Time (Myr) 

Figure 12. The change in total magnetic energy of all the runs. 
The magnetic energy is calculated only considering cells whose den- 
sity is greater than 1.01 times the initial density. The dashed curve 
represents the mhd curve divided in half. 

fast magnetosonic wave leaves the grid, even though in 
reality magnetic energy is constantly being injected into 
the GMC. 

In order to account for only the magnetic energy in- 
jected into the GMC by the HII region and not the initial 
magnetic energy present in each cell, we calculate the 
total change in magnetic energy, AEb, where AEb = 
{{B,f + {Byf + {B^f)/2-{{B,,f + {By„f + {B^„r)/2, 
where , By^ , Bzq is the initial magnetic field strength 
in the x, y, and z directions, respectively. We find that 
this quantity is very similar for all the blister runs. Fig. 
13 shows the specific change in magnetic energy AEb /M. 
The curves all decline with time at late times, which im- 
plies that the rate at which the magnetic energy is in- 
jected is slower than the rate at which mass is swept up 
by the ionization front for all the runs. The symmetric 
HII regions have both more total magnetic energy and 
specific magnetic energy than their blister counterparts, 
although at late times the difference in the specific mag- 
netic energy is very small. 
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Figure 13. The change in specific magnetic energy. 
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Figure 14. The total energy of ah the runs. The dashed curves 
represent the hydro and MHD runs divided by a factor of two. 



3.3.4. Total Energy 

In Fig. [T4j we plot the total energy, that is, kinetic 
plus magnetic energy, for all the runs. For the hydro and 
blister- hydro runs, there is no magnetic energy of course, 
so this just consists entirely of the kinetic energy. It is 
clear that magnetic effects are of great importance with 
respect to the total energy of the cloud. By 5 Myr, the 
blister-mhd run has about twice as much energy as the 
blister-hydro run, and the blister-mhd-vert and blister- 
mhd-45 runs have about 30 % more energy. We plot the 
specific total energy in Fig. [15] The specific total energy 
depends most strongly on the magnetic field orientation 
and not on whether it is a blister or symmetric case. 
There is no advantage in going from hydro to mhd in 
the symmetric case, but there is an advantage in doing 
so for the blister case at least for some magnetic field 
orientations. 

It is instructive to compare the kinetic energy lost to 
the magnetic energy gained in going from blister-hydro 
to blister-mhd. The blister-mhd case has 6 x 10^^ erg 



less kinetic energy at 5 Myr, but it also has ^ 3 x 10^^ erg 
more magnetic energy. Thus the gain in magnetic energy 
outweighs the loss of kinetic energy by a factor of 5. As 
a result, we find that in general that magnetized HII 
regions deliver significantly more energy to dense clouds 
than their pure hydrodynamic counterparts. 

4. DISCUSSION 

Here we investigate the origin of the differences be- 
tween the analytic approximations and the numerical re- 
sults that we described in the previous section. 

The validity of the approximations used to derive the 
analytic approximation can be checked through a simple 
calculation. Using Eqn. ([lO| to estimate the HII region 
intern al d ensity, and a value for the ionizing luminosity 
from §2.3[ we find that the ratio of ambient pressure to 
HII region pressure is ~ 4.4% at 1 Myr and ~ 12.3% 
at 3 Myr. Thus although the approximation of HII re- 
gion pressure dominance works quite well at early times, 
by 3 Myr it is not quite as accurate, and this discrep- 
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Figure 15. The specific total energy of all the runs. 



ancy can cause some of the observed flattening of the 
bUster-hydro and spherical-hydro curves late in the run, 
although clearly not enough to account for the substa- 
tional flattening in the blister-hydro run observed in Fig. 

m 

In reality this flattening is probably due to the defi- 
nition used to calculate the radius of the blister-hydro 
curve. The radius is a bit misleading because it turns 
out that the spherical portion of the blister-hydro shell 
travels significantly slower than the analytic solution pre- 
dicts. This is shown in Fig. [TBj where we have calculated 
the radius of the blister-hydro shell by only considering 
gas along the y = z = line. The spherical part of 
the shell expands just a bit faster than the hydro shell. 
This implies that the slivers of the shell that extend in 
the y direction along the interface (as seen in Fig. l4| 
contribute significantly to the radius of the blister hydro 
curve. In the early to mid-stages of the simulation the 
slivers expand faster in the y-direction than the spheri- 
cal part of the shell expands in the radial direction, but 
towards the later stages the expansion of the slivers in 
the y-direction slows down dramatically and the curve 
becomes flatter than the analytic solution. 

While this may explain the shape of the curve, it does 
not explain why the blister-hydro shell seems to expand 
at a rate comparable to the hydro shell, rather than 
roughly 20% faster as predicted by the analytic approx- 
imation. One possible explanation for this discrepancy 
is that the analytic solution for the blister case assumes 
that the density is the same as in the symmetric HII re- 
gion. The D-type ionization front travels below the sound 
speed of the HII region, so gas inside the symmetric HII 
region has time to spread out and achieve a uniform den- 
sity. This is not quite true for the blister case since the 
gas is free to escape from the HII region, so that there is 
a non-uniform density distribution inside the blister-type 
HII region. We check this by comparing the density just 
inside the shell for both the blister-hydro and hydro runs 
(Fig. 17). The blister-hydro density is about 15% less 
than tne hydro density for the entire run. If we multiply 
pii in Eqn. |Al| by a factor of 0.85 in the blister case to 
account for this effect, we find the predicted difference in 
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Figure 16. The radius for the blister- hydro run computed by only 
considering gas along the y = z = line. The blue-dashed curves 
are the same as in Fig. [s] 
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Figure 17. The density just inside the dense shell of the hydro 
and blister-hydro HII regions plotted over time. 

shell radius and speed between the symmetric and blis- 
ter cases drops from 20% to 15%. Thus the incorrect 
assumption of a uniform density in the blister case ac- 
counts for about a quarter of the discrepancy. The rest 
is likely due to a failure of the assumption of hemispher- 
ical symmetry. Comparing Fig.'s |8] and [16) we see that 
the mass-averaged radius considering all angles is ~ 50% 
larger than the value along the y = z = line even at 
very early times, before significant tails form. Thus the 
shell is only very roughly hemispherical. These results 
show that the blister HII regions are much more complex 
than their symmetric counterparts, and that the analytic 
solution for the blister case cannot use the same simpli- 
fying assumptions that work quite well for the symmetric 
case. 

5. CONCLUSIONS 

We have performed the first numerical study of a blis- 
ter type HII region expanding into a magnetized medium. 
We draw the following conclusions: 
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(i) Although the kinetic energy of the magnetized runs 
is lower than that of their hydrodynamic counterparts, 
they have much more total energy since the kinetic en- 
ergy lost in going from hydro to MH D is many times 
less than the magnetic energy gained ( ^3.3.4 ), and hence 
could possibly be more efficient at drivmg turbulence. 
It is not entirely clear how efficient the injected mag- 
netic energy is at driving turbulence compared to the 
kinetic energy, but studies of Alfven waves decaying into 
turbulence in various astrophysical environments (from 
GMCs to the solar wind) suggest that the process is 
likely to be very efficient. A circular Alfven wave devel- 
ops a "decay" instab ility which ultimately leads to decay 
into turbulence (see McKee & Ostriker (2007) and ref- 
erences therein). This decay into turbulence requires an 
initial directional imbalance in the Alfven waves, which 
we have in our simulations since the waves are all left- 
propogating. It is likely that the magnetic energy added 
to the cloud by including the effects of MHD is at least 
as important, if not more important, than the kinetic 
energy for any type of HII region. Therefore it is impor- 
tant to include the effects of MHD in future studies of 
star formation. 



(ii) A blister-type HII region expands into a cloud 
faster than the corresponding symmetric case, but not 
by as much as predicted by some analytic approxima- 
tions. 

(iii) The total energy is greatest in the symmetric case, 
so HII regions of this type make the greatest contribu- 
tion to the total energy budget of a cloud. However, since 
GMCs are turbulent and have a filamentary morphology, 
most new born stars are likely to be near the edge, so the 
blister scenario should be more common than the embed- 
ded one. Nonetheless, our simulations confirm that, in 
the presence of a magnetic field, even blister-type HII 
regions can inject significant energy into the dense parts 
of molecular clouds. 
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APPENDIX 

DERIVATION OF THE SHELL EXPANSION RATE 

Here we derive the analytic solution of the radius of the shell in the spherical and blister cases. We assume ionization 
balance, consider the density pn constant inside the HII region, and use momentum conservation: ^ = J^, where 
V is the momentum of the shell and T is the force applied to it by the matter inside the HII region. The mass of 
the shell is Mgh = (4, 2)7rr^po/3 since most of the mass inside the Stromgren Sphere is contained in the shell to good 
approximation, where 4 and 2 are the coefficients for the symmetric and blister cases, respectively. Using ionization 
balance and Eqn. (p^ol) we can write the density inside the HII region 



Pii 



V47ra(^) 



1/2 



-3/2 



(Al) 



1/2 



-3/2 



(A2) 



and hence the pressure (thermal plus ram) inside the HII region 

P = il,2)puch = il,2)ch 

where the coefficient of 2 for the blister case represen ts the ram pressure of material rocketing off the inside of the 
dense shell as it is ionized ( Krumholz fc Matzner|20Q9 ). Now we use momentum conservation to arrive at the equation 
of motion: 



[(4, 2)7rrVor/3] =:F = PA = Anr'ch (^|y 



dV _ d 
dt dt 

where A is the surface area of the sheh. 
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This ODE admits a similarity solution of the form r (xV ^ and with some algebra one can show that 



and 



where tg = Vsjcu. 
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(spherical). 
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(blister). 



(A3) 



(A4) 



(A5) 



(A6) 



REFERENCES 



14 



Gendelev & Krumholz 



Arthur, S., Henney, W., Mellema, G., Colle, F., &; Vazquez-Semadeni, E. 2011, MNRAS, 541 
Crutcher, R. M. 1999, ApJ, 178 

Dale, J. E., Bonnell, 1. A., & Whitworth, A. 2007a, MNRAS, 375, 1291 
Dale, J. E., Clark, P. C, & Bonnell, I. A. 2007b, MNRAS, 377, 535 
Elmegreen, B. & Lada, C. 1977, ApJ, 214, 725 

Fukui, Y., Kawamura, A., Minamidani, T., & et al. 2009, ApJ, 705, 144 

Goldbaum, N., Krumholz, M., Matzner, C, & McKee, C. 2011, ApJ, 738, 101 

Gritschneder, M., Naab, T., Walch, S., Burkert, A., & Heitsch, F. 2009, ApJ, 694 

Henney, W., Arthur, S., Colle, F., & Mellema, G. 2009, MNRAS, 398, 157 

Israel, F. 1978, A&A, 70, 769 

Kahn, F. 1954, BAN, 12, 187 

Koyama, H. & Inutsuka, S. 2002, ApJ, 564 

Krumholz, M., Matzner, C, & McKee, C. 2006, ApJ, 653, 361 

Krumholz, M. & Matzner, D. 2009, ApJ, 703, 1352 

Krumholz, M., Stone, J., & Gardiner, T. 2007, ApJ, 671 

Krumholz, M. &; Tan, J. 2007, ApJ, 654, 304 

Mackey, J. & Lim, A. 2010, MNRAS, 412, 2079 

Matzner, C. 2002, ApJ, 566, 302 

McKee, C. &; Ostriker, E. 2007, ARA&A, 45, 565 

McKee, C. &; Williams, J. 1997, ApJ, 476, 144 

Mellema, G., Arthur, S., Henney, W. J., I.T., I., & Shapiro, P. R. 2006, ApJ, 647, 397 

Osterbrock, D. 1999, Astrophysics of Gaseous Nebulae and Active Galactic Nuclei (Mill Valley: University Science Books) 

Pellegrini, E., Baldwin, J., Brogan, C., Hanson, M., Abel, N., Ferland, G., Nemala, H., Shaw, G., & Troland, T. 2007, ApJ, 658, 1119 

Spitzer, L. 1978, Physical Processes in the Interstellar Medium (New York: Wiley) 

Stone, J., Gardiner, T., Teuben, P., Hawley, J., & Simon, J. 2008, ApJS, 178 

Stromgren, B. 1939, ApJ, 89, 526 

Tenorio-Tagle, G. 1978, A&A, 71, 59 

Whitworth, A. 1978, MNRAS, 186, 59 

Whitworth, A., Bhattal, A., Chapman, S., Disney, M., & Turner, J. 1994, A&A, 290, 421 
Williams, J. & McKee, C. 1997, ApJ, 476, 166 



